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^\ , We introduce a microscopic model of a lipid with a charged headgroup and flexible hydropho- 

0^ ' bic tails, a neutral solvent, and counter ions. Short-ranged interactions between hydrophilic and 

0^ I hydrophobic moieties are included as are the Coulomb interactions between charges. Further, we 

include a short-ranged interaction between charges and neutral solvent, which mimics the short- 
f*"! ' ranged, thermally averaged interaction between charges and water dipoles. We show that the model 

^ , of the uncharged hpid displays the usual lyotropic phases as a function of the relative volume fraction 

of the headgroup. Choosing model parameters appropriate to dioleoylphosphatidylethanolamine in 
^SJ ' water, we obtain phase behavior which agrees well with experiment. Finally we choose a solvent 

concentration and temperature at which the uncharged lipid exhibits an inverted hexagonal phase 

and turn on the headgroup charge. The lipid system makes a transition from the inverted hexago- 
^^ I nal to the lamellar phase which is related to the increased waters of hydration correlated with the 

increased headgroup charge via the charge-solvent interaction. The polymorphism displayed upon 

variation of pH mimics that of the behavior of phosphatidylserine. 
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I. INTRODUCTION 



Biological lipids in solution display several different lyotropic phases, and the implications this may have for bio- 
logical function has been a subject of speculation for many years (Cullis et al., 1985; de Kruijff 1997). Lipid phase 
behavior depends upon several factors, some of which are intrinsic to the lipid architecture itself. For example, an 
increase in the length of the hydrocarbon tails brings about transitions from lamellar. La, to inverted hexagonal, Hn, 
phases (Seddon, 1990), while an increase in the volume of the headgroup brings about the reverse (Gruner, 1989). 
Other factors regulating phase behavior are externally controlled, such as temperature, solvent concentration, and 

^1- ■ solvent pH (Hope and Cullis, 1980, Seddon et al., 1983, Bezrukov et al., 1999). It is these factors which arc the focus 

^^ ', of this paper. 

OQ ' Lipid phase behavior has been addressed extensively by the construction of phenomenological free energy functions 

^^ which contain terms describing, inter alia, bending, hydration, and interstitial energies (Helfrich, 1973; Kirk et al., 
1984; Rand and Parsegian, 1989; Kozlov et al., 1994). Such approaches, which obtain their several parameters from 
experimental measurement of various quantities, are quite useful, particularly in correlating phase behavior with 
other thermodynamic properties. Nonetheless, it would clearly be desirable to derive all thermodynamic quantities, 
including the phase behavior, by applying statistical mechanics to a microscopic model of the system. In addition to 
simplifying the description considerably, such approaches would correlate phase behavior with the the architectural 

'^ ■ properties of the lipid itself and its solvent. 
^ ' Analytic, mean-field, approaches of statistical mechanics have been applied to anhydrous lipids to investigate 

O ' behavior of increasing complexity. Such methods have been combined with realistic models of lipid tails to determine 
how the hydrocarbon chains pack in aggregates and in bilayers (Marcelja, 1974; Gruen, 1981 and 1985; Ben-Shaul 
et al. 1985; Fattal and Ben-Shaul, 1994). Results for the bilayer are in good agreement with molecular dynamic 
simulation (Heller et al., 1993). These methods have shown that in a neutral, anhydrous system, the entropy of the 
lipid tails always favors the Hu over the La phase, and that a change in area per headgroup could bring about a 
J^ transition between them (Steenhuizen et al., 1991). 

Aggregates, such as the lipid bilayer, in the presence of solvent have also been considered within the mean-field 
approach applied to lattice models (Leermakers and Scheutjens, 1988). In addition to the tails, one must now model 
the solvent and the headgroups, and phosphatidylcholine and phosphatidylserine headgroups are among those which 
have been described (Meijer et al., 1994). The method is flexible and has been applied to many different systems, 
including bilayers with trans-membrane guest molecules (Leermakers et al., 1990). Results are quite good, with the 
exception that the local volume fraction of solvent inside the bilayer is rather large, several orders of magnitude 
greater than that observed in experiment (Jacobs and White, 1989). Lattice models, however, are not well-suited to 
the description of transitions between phases of different symmetry. 

Recently similar methods were applied to a system of solvent and monoacyl lipid embedded in a continuous space, 
and the phase diagram was obtained by solving the mean-field theory exactly (Miiller and Schick, 1998). It displayed 



both La and Hjj phases, so that the transition between them could be studied as a function of hpid architecture. 
The dependence of the transition on the architectural parameters, length of tail and volume of headgroup, was that 
observed in experiment. However, the fraction of solvent within the bilayers was again too large. 

In this paper, we introduce a computationally more tractable model of a lipid than that employed by Miiller and 
Schick, one whose hydrocarbon tails arc modeled as flexible chains rather than within the rotational isomeric states 
framework (Flory, 1969; Mattice and Sutcr, 1994) employed earlier. We first study the model with an uncharged head- 
group. Its phase behavior, both with respect to variations in architecture and variations in solvent concentration, is as 
expected, and in agreement with experiment. In particular, choosing model parameters appropriate to dioleoylphos- 
phatidylethanolamine (DOPE), we obtain a phase diagram similar to that observed (Gawrisch et al. 1992, Kozlov et 
al. 1994) . We extract the variation with temperature and solvent concentration of the lattice parameter of the inverted 
hexagonal phase, and compare it to experiment (Tate and Gruner, 1989, Rand and Fuller, 1994). The agreement is 
excellent. We also find that the concentration of solvent within the bilayer is vanishingly small. We then allow the 
headgroup to be negatively charged, introduce counter ions into the system, and include the Coulomb interaction 
between all charges. We show that the Coulomb repulsion between the headgroups does not tend to stabilize the La 
phase at the expense of the Hn, but has the opposite effect. Lastly we turn on a short-ranged interaction between 
charges and neutral solvent, an interaction which models the thermally averaged interaction between charges and the 
dipole of water. We then find that as the charge on the headgroup is turned on, the La phase is stabilized with 
respect to the Hu. In effect, as the charge on the headgroup increases, so too do the waters of hydration. In addition, 
the counter ions which are attracted to the headgroup are also enlarged by their own waters of hydration. It is the 
totality of these waters which effectively increases the headgroup volume and therefore stabilizes the lamellar phase. 

The paper is organized as follows. In the next section, we introduce the model for the charged lipid, the solvent, 
and counter ions, specify all the interactions between them and set up the partition function of the system. In section 

III, we derive the self-consistent field theory for it. At the heart of the theory are four self-consistent equations for the 
electrostatic potential of the system and the three effective fields which determine the headgroup, tail, and solvent 
densities. One of these self-consistent conditions is simply the non-linear Poisson-Boltzmann equation. In section 

IV, we expand all functions of position into a complete set of functions having a specified space-group symmetry, 
and rewrite the self-consistent equations in terms of the coefficients of these expansions. These equations are solved 
numerically, and the free energies of the various phases computed. A comparison of the free energies yields the 
phase diagram. In Section V, we present the phase diagram for the neutral lipid as a function of temperature and 
one architectural parameter. We include here only the classical phases, lamellar, inverted and normal hexagonal, 
and inverted and normal body-centered-cubic, as well as the disordered phase. For the remainder of the section, we 
choose an architecture such that the anhydrous, neutral lipid orders into the Hjj phase. Results for the system in 
the presence of a neutral solvent, along with comparisons to experiment, are presented next. 

In Section VI, we consider the charged lipid. We choose a water concentration such that the neutral lipid remains 
in the Hu phase. By varying the counter ion concentration, we turn on the charge on the headgroup, and thus all 
Coulomb interactions, but we keep the short-ranged interaction between charges and neutral solvent set to zero. We 
show that the Hu phase is stabilized with respect to the La rather than destabilized, a result completely analogous to 
that of polyclectrolytes (Nyrkova et al., 1994). Lastly we present our results for the case in which, in addition to the 
Coulomb interaction, a short-ranged interaction between charges and solvent is turned on. Here we find that the La 
is indeed stabilized with respect to the Hjj phase, in agreement with experiment (Hope and CuUis, 1980, Bezrukov 
et al., 1999). 

II. THE MODEL 

We consider a system composed of charged lipids, neutral solvent, and counter ions in a volume V. There are ul 
lipids, each of which consists of a head, with volume Vh, and two equal length, completely fiexible, tails each consisting 
of N segments of volume Vt- Each lipid tail is characterized by a radius of gyration Rg = (A^a^/6)^/^, with a the 
statistical segment length. The heads carry a negative charge —eQh- The solvent consists of Us neutral particles of 
volume Vs, while the ric counter ions have charge +e and negligible volume, Vc = 0. There are five dimensionless 
densities which totally specify the state of the system; the number density of the headgroups, <I>^, of the tail segments, 
$4, and of the solvent, $s, and the charge density of the headgroups, ePh, and of the counter ions, ePc- They can be 
written as 

<l;,(r)=w^^(5(r-r,(l/2)), (1) 

1=1 



"1 „1 
<|.,(r)=t;;,V / 6{r - ri{s))ds, (2) 



<i.,(r)-r;;,£<5(r-R,,,), (3) 

A(r)--«,.^0M'^(r-r,(l/2)), (4) 

ric 

A(r)-i;„;^<5(r-R,,,). (5) 

1=1 

We have chosen v^ as a convenient volume to make aU densities dimensionless. In the above, R^.j is the position 
of the j'th solvent particle, and RcA the position of the i'th counter ion. The configuration of the Tth lipid is 
described by a space curve r;(s), where s ranges from 0, at the end of one tail, through s = 1/2 at which the head 
is located, to s = 1, the end of the other tail. The nominal probability that the charge on the headgroup of the Z'th 
lipid, — e Qh,h is equal to — e or is p or 1 — p respectively. As we model the case in which charges can associate or 
disassociate from the headgroup, it will be necessary to average the partition function of the system with respect to the 
charge distribution. This corresponds to an annealed distribution in the nomenclature of Borukhov (Borukhov ct al., 
1998). The concentrations of lipid, solvent, and free counter ions are controlled by chemical potentials. In particular, 
increasing the number of free, positive, counter ions implies, by charge neutrality, an increase in the negative charge 
on the headgroups, and thus corresponds to an increase in the pH of the system. 

The interactions among these elements are as follows. First there is a repulsive, contact interaction between 
headgroup and tail segments, and also between solvent and tail segments. The strength of the interaction is kTvhXi 
where k is Boltzmann's constant and T the absolute temperature. Second there is the Coulomb interaction between 
all charges. The dielectric constant of the solvent is denoted e. Finally there is a contact interaction between all 
charges and the neutral solvent whose strength is kTviiX. This is to model the short-ranged, thermally averaged, 
interaction between charges and the dipole of water, an attractive interaction which decreases like r^^ and is of 
strength e^v? /Qe^kT , where u is the dipole moment of water (Israelachvili, 1985). Thus the energy per unit volume 
of the system, E/V , can be written 

^f [*;„*„*.. ft. ftl = 2xA'y*|4k(r) + 4.(r)]4,(r) 

f dr ~ 

^A /-$,(r)[P,(r)-P,(r)], (6) 



where 



V 






VhtkT 

is a dimensionless measure of the strength of the Coulomb interaction. The grand partition function (Matsen, 1995) 
of the system is 



ni ric ns r ' ^ ^ 

^ ^ift^! / n T'^^T'Q^^^ n '^i^- n ^^s, 

ni,nc,ns l — l i—1 j — 1 

exp {-E[^h, 4t, $„ Ph, Pc]/fcr} 6{l -^h- 7s$. - lA). (8) 



Here jT^ri denotes a functional integral over the possible configurations of the Tth lipid and in which, in addition to 
the Boltzmann weight, the path is weighted by the factor 7'[rj ;(s); 0, 1], with 



V[r,si,S2] =:7Vcxp 



SRI Js. '' ds I . 



(9) 



with Af an unimportant normalization constant. The notation J'DQh^i denotes an integral over the probability 
distribution of the charge on the headgroup of the Tth lipid. We have enforced an incompressibility constraint on 



the system with the aid of the delta function S{1 — <bh — 7s^s — Jt^t), where 7^ = Vg/vh, and jt = '^Nvt/vh- The 
latter parameter is the lipid architectural parameter. The relative volume of the headgroup with respect to that of 
the entire molecule is 1/(1 + 7t)- 

The model is now completely defined. The solvent is specified by 7s, its volume per particle relative to that of 
the headgroup, and the architecture of the lipid is characterized by 74. There are three interactions, hydrophobic- 
hydrophilic, charge-charge, and charge-solvent, whose strengths are given by x, /?*, and A respectively. The external 
parameters are the temperature, conveniently specified in terms of a dimensionless temperature T* = (2x-/V)~^, the 
fugacity of the solvent, z^, and the fugacity of the free counter ions, Zc, which, by charge neutrality, controls the 
charge on the lipid headgroups. The characteristic length in the system is the radius of gyration, Rg. In the next two 
sections we derive the self consistent field theory for the model, first in real space, and then in the following section, 
in Fourier space. 

III. THEORY: REAL SPACE 

Evaluation of the partition function of Eq. H is difficult because the interactions are products of densities, each 
of which depends on the specific coordinates of one of the elements of the system. This dependence is eliminated 
in a standard way. We illustrate it on $/i(r) which, from its definition in Eq. H, depends on the coordinates of the 
headgroup, r;(l/2). One introduces into the partition function the identity 

1= fv^hS{^h-^h), 
= J V<^hT^W,,expij- J Wi,{r)[<i>h{r) - ^h{r)]dr\ , (10) 

in which <&/i(r) does not depend on any specific coordinates of one of the elements of the system, but is simply a 
function of r. The integration on Wh extends up the imaginary axis. Inserting such identities for the five densities 
<l>/j, <l>t, $s, Pft, and Pc, and a similar identity for the delta function expressing the incompressibility condition, one 
rewrites the partition function, Eq. 0, as 



Z = / V^h'DW,,V<Pt'DWtV^s'DWsVPhVUhVP,VU,VE 

X eMziQdWh, Wu Uh] + z,Q,[U,] -f z^QsWs] - E[<^h, <^u *., Ph,Pc]/kT} 

X exp |i- /" [Wh-^h + Wt-^t + Ws<^s + UhPh + U,Pc + S(l - $^ - 7,$, - 7t$t)] drX 



(11) 



where 



Qi[Wn,Wt,Uh]^ j'Dvi'bQhexv[-Wh{vi{l/2)) + QnUn{vi{l/2))~ j dsWt{vi{s))\ , (12) 

is the partition function of a single lipid in external fields Wh, Wt, and Uh, 

Qc[Uc] = f dRceM-UcC^c)], (13) 

is the partition function of a single counter ion of unit positive charge in an external potential Uc, and 

Q.[WA = /"dR,exp[-I¥,(Rs)], (14) 

is the partition function of a single solvent particle in the external field Wg. It is convenient to shift the zero of all 
chemical potentials so that zi — > l/v^, Zc -^ Zc/vh, and Zg — > Zg/vh- The partition function, Eq. |ll|, can then be 
written in the form 

Z= f V^hVWhV^tVWtV^gVWgVPnVUhVPcVUcV^ exp[-f]/fcr] , (15) 

with 



... ^^ ^ _ a[m..»'.t/,.i _ ^^^ap _ „M!iJ + j;^m*,.,*.,*„p,..pj 



kTV V "■ V V kTV 

f dr 
- J yiWh^h + Wt^t + Ws<i>s + UhPh + U,P, + S(l - $„ - 7,$, - 7t$,)]- (16) 

No approximations have been made to this point. What has been accomphshed is a rewriting of the partition function 
from a form, EqH, in which all entities interact directly with one another, to a form, Eqs. 15 and Hq, in which they 
interact indirectly with one another via fluctuating fields. Although the integrals in Eq. [I^ over $^, $(, $s, P/i, Pc 
and S could all be carried out, as they are no worse than Gaussian, the integrals over the fields Win Wt, Wg, Uh, 
and Uc cannot. Therefore we employ the self-consistent field theory in which we replace the integral in Eq. g^ by 
its integrand evaluated at its extremum. The values of Wu, <i>;i, etc. which satisfy the extremum conditions will 
be denoted by the corresponding lower case letters Wh, and (ph, etc. The equations which determine them are six 
self-consistent equations for the six fields Wh-, Wt, Wg, uu, Uc, and ^. They are 

Wh{v) = 2xNM-r)+(.{v) (17) 

wt{v) = 2xNiMr) + Mr)) + 7*^(0 (18) 

Wgir) = 2xNMr) - A(pc(r) - p„(r)) + 7,e(r) (19) 

Uhir) +Uc{r) 
w(r) = 2 

^(r_ fdv' ph{v') + Pc{v') , . 

in J Rl |r-r'| ^ ' 

^XMr) (21) 

l = 0h(r)-)-7t0t(r)-)-7,0,(r). (22) 

As the field ^ is easily eliminated, the six equations readily reduce to five. The simplicity of Eq. G^ reduces this, in 
practice, to a set of four equations. The five densities (ph-, 0t, P's, Ph, and pc are functionals of all of the above fields 
except ^, and therefore close the cycle of self-consistent equations: 

(j>h{r)[wh,wt,Uh] = J — ^^ 23) 

, , xr 1 SQl[Wh,Wt,Uh] ,„.. 

<j>t{r)[wh,wt,Uh] = ^^ ., (24) 

(f>s{r)[ws] = ~Zs- — -— (25 

ows[r) 

= Zsexp[-Ws{r)] (26) 

/NT 1 SQi[wh,Wt,Uh] , . 

Ph{r)[wh,wt,Uh] = J — 1-^ (27) 

Pc r M = - 28 

duc(r) 

= Zcexp[-Mc(r)] (29) 

The density 0/i(r) is simply the expectation value of <I'/i(i") in the single lipid ensemble. Similar interpretations follow 
for the other densities. Note that one of the self-consistent equations, Eq. ISO, is simply the non-linear Poisson- 
Boltzmann equation, and u(r) the electric potential. 

With the aid of the above equations, the mean- field free energy, f^m/, which is the free energy function of Eq. 11^ 
evaluated at the mean- field values of the densities and fields, can be put in the form 

kT 

- Q,nf = {Ql[wh,Wt,Uh] + ZcQc[Uc] + ZsQs[Ws]) + E[(j)h, P't, (t^s, Ph, Pc], (30) 

Vh 

= kT{ni + nc + Ug) + E[(j)h,4>t, (t>s, Ph, Pc], (31) 

with E given by Eq. pi The thermodynamic potential, O, is that appropriate to an incompressible system calculated 
in the grand ensemble; the negative of the osmotic pressure multiplied by the volume. Thus the above equation states 



that the osmotic pressure is the sum of the ideal, partial osmotic pressures plus a correction due to the interactions. 
Within mean field theory, this correction is simply the energy per unit volume of the system. 

We now specify that the charges in the system can associate with or disassociate from the headgroup in response 
to the local electrostatic potential. This implies that the partition function of a single lipid, Qi is to be averaged over 
the nominal charge distribution that Qh — 1 with probability p, and Qh — with probability I ~ p (Borukhov et al., 
1998). The consequence of this averaging is that Q;[w;i, Wt,u/i] of Eq. [l^ becomes 

Qi[wh.,eff,wt]^ i>T^ie^P<~Wh,effiri{l/2))- dswt{ri{s)) [ , (32) 

where 

Wh,effir:) = Whir) - In / VQh exp[QhUh(r)] (33) 

= Whir) - ln[l +p(expK(r)] - 1)]. (34) 

Although this appears to introduce an unknown parameter p into the problem, the condition of charge neutrality, 

J dr[phir) + pcir)] ^ 0, (35) 

relates this parameter to the fugacity of the counter ions, Zc- In practice, we use this fugacity to control the pH and 
the amount of charge on the lipids. 

There remains only to specify how the single-lipid partition function is obtained. One defines the end-segment 
distribution function 

Wh,eff{riit))6it-l/2)+wtiriit)))}, (36) 



qir, s) ^ Vri{s)d{r - r,(s)) exp <^ - 


r ( 

Jo V 


1 .dRit) 2 

.^Rg dt ' _ 


which satisfies the equation 




^'^(^'^)=2 7?2X7 


"^nir. ,s^ - 


Uiii, „ff(r)()(s ^ 



2)+Wtir)]qir,.s), (37) 

with initial condition 

'?(r,0) = l. (38) 

The partition function of the lipid is then 

Qi^ fdrqir,!). (39) 

From this expression for the single-lipid partition function and Eqs. p3|, p4|, and p^ one obtains expressions for the 
local density of the lipid heads 

(/),,(r) = e^p[-Wh.effir)]qir, ^-)qir, --), (40) 

of the lipid tails 

(j)tir) = / dsg(r,s)q(r,l-s) (41) 

Jo 

and of the charge density on the lipid heads 

f,.\ pexp[uhir)] 

Phir) = -— — -, -, — -r--. -T0h r). 42) 

l-|-p(exp[uft(r)J - 1) 

To summarize: there are four self-consistent equations to be solved for the fields Wh, wt, Wg, and electrostatic 
potential u. These equations, obtained from simple algebraic manipulation of Eqs. n^ to E2l can be taken to be 



jtwhir) - wtir) = 2xN[jtMr) ~ Mr) - Mr)], (43) 

jsWhir) - ws{r) = 2xN{js-l)Mr) + KPc{r)'Ph{r)), (44) 

1 = Mr)+ltMr)+lsMr), (45) 

RlvMr)^-(3*{ph{r)+Pc{r)). (46) 

Note that we have chosen here to write the Poisson-Boltzmann equation, Eq. EOl in its local, rather than its integral 
form. When the four fields are known, the corresponding densities follow from Eqs. |2^, |2^, ^, |4|, and ^. 

Rather than attempt to solve these equations in real space, a difficult task for the periodic phases in which we are 
interested such as iJ//, we recast the equations in a form which makes straightforward their solution for a phase of 
arbitrary space-group symmetry (Matsen and Schick, 1994). 

IV. THEORY: FOURIER SPACE 

We note that the fields, densities, and the end point distribution function depend only on one coordinate r. Therefore 
in an ordered phase, these functions reflect the space-group symmetry of that phase. To make this symmetry manifest 
in the solution, we expand all functions of position in a complete, orthonormal, set of functions, /i(r), i = 1,2,3..., 
each of which have the desired space group symmetry; e.g. 

Mr)^Y.^hMr), (47) 



1 

V 



M = T7 / dvh{v)f,{v). (48) 



Furthermore we choose the fi (r) to be eigenfunctions of the Laplacian 






VV.(r) = -7^/'(r), (49) 



where _D is a length scale for the phase. The functions for the lamellar phase are clear. They can be taken to be 

/i(r) = 1, (50) 

f^{v) = \/2cos[27r(i - l)x/D], i > 2, (51) 

Expressions for the unnormalized basis functions for other space-group symmetries can be found in X-ray tables 
(Henry and Lonsdale, 1969) as they are intimately related to the Bragg peaks. In the tables cited, those for the 
hexagonal phase, space group (p6m) can be found on page 372, and that of the bcc phase, space group (ImSm) on 
page 524. 

The four self-consistent equations become 

"ftWh;i - wt-i = 2xNYytM - Ml - (I^s-a], (52) 

"fsWh-i - Ws-i = 2xAf(7s - l)0t;i -I- X{pc-i - ph-i), (53) 

Sl,i = (t>h:i + 7t0t;'i + ls4>8-i, (54) 

\ 7-)2 
-Jyf'^^ = P*iPh;i + Pc-i) (55) 

To obtain the partition functions and densities, we proceed as follows. For any function G'(r), we can define a 
symmetric matrix 

{G)^j^^ j h{v)G{v)f,{v)dv (56) 

Note that {G)ii = {G)ii ~ Gi, the coefficient of /i(r) in the expansion of G{r). Matrices corresponding to functions 
of G'(r), such as 



ie%, - 



^^j h{v)e^^^^f,{v)d{v), (57) 



are evaluated by making an orthogonal transformation which diagonalizes {G)ij. With this definition, Eqs. |2^ and 
|29| yield the solvent density and counter ion charge density 

q^s:,=Zs{e~^')^,, (58) 

Pc■^ = Zc{e-''^)^^^, 

= Zc(e-("-"=))^^. (59) 

(60) 
To obtain the remaining densities, we need the end-point distribution function. From Eq. pj we obtain 

^ ^-J2lA., + {wn,eff\A^ - l/2)]g,(s), (61) 

3 

2Rl 
^-1 = -^^-^-0 + {wth. (62) 

with initial condition gi(0) = ^^^i. The solution of this equation is 

g,(,s) = (e-^0, 1 , if s < 1/2 

-m^~-''-").,{--^")^,^ ^-1/2 

= E {^-^^'-'''') . . (e-^"'-»),, {e-^'^) , -^ > 1/2. (63) 

From this, the remaining densities follow from Eqs. Efl, O, and 42: 



jkl 



J2{e--^^^ff)^^T,,Ml-)<i4-), (64) 



., = / ds V'ry7,qj(s)(7A;(l - s), 



(65) 



pgUh 



■J 
with 



PM = -E ! + „(,«. _l) )^./''u- (66) 



The mean-field free energy, Eq. pU takes the form 



r.^fe^-^ //.(r)/,(r)A.(r). (67) 



with the mean-field energy being given by 

kTV 



kTV 

^mf = (0t;l+0s;l+Pc;l)+£^, (68) 



E^ 



Vh 



^[2xN{(ph-i + (l)s-s)4>t:i + 7^{ph-i + Pc-i)ui - \{pc-i - Ph-i)4>s;i]- (69) 



We have expressed the Coulomb energy as a product of the charge densities and electrostatic potential. Note that this 
free energy still depends parametrically on D, the length scale of the phase, so that the value of D which minimizes 
it must be determined. Once this is done, we compare the free energies obtained for phases of different space-group 
symmetry, and thereby determine the phase diagram of our model lipid system. 

The infinite set of self-consistent equations, Eqs. |5^ to |5^ must be truncated in order to be solved numerically. We 
have employed up to 50 basis functions. This truncation is sufficient to ensure, for T* > 0.03 and 1/(1 -I- 7*) < 0.66, 
an accuracy of 10"'* in the free energy Vh^mf/kTV. As noted, one must also determine the length scale which 



minimizes the free energy. This is usually straightforward as there is a single well-defined minimum for a phase of 
given symmetry at given thermodynamic parameters; temperature, and chemical potentials. Were there more than 
one minimum, this would reflect a tendency for the system to phase separate into two phases with the same non-trivial 
space-group symmetry, an extremely unusual occurrence. The single minimum that one finds normally is sharp, that 
is, the free energy varies rather rapidly with D. Only in cases in which phases are greatly swollen is the minimum 
extremely shallow and difhcult to locate. 

V. RESULTS: THE NEUTRAL LIPID 

We first apply our method to a neutral lipid. We show here the phase behavior of the neutral lipid, in the absence 
and in the presence of solvent. Figure 1 shows the phase diagram of the pure lipid as a function of the dimensionless 
temperature T*, and the architecture of the lipid. The latter is characterized by the single parameter 1/(1 -1-74) which 
is the relative volume of the headgroup to that of the entire lipid. It is analogous to, but not the same as, the single 
parameter used by Israelachvili to characterize the geometry of lipids (Israelachvili, 1985). Shown are the lamellar 
phase. La, the normal and inverted hexagonal phases, Hj and Hn, and the normal and inverted body-centered cubic 
phases bcci and bccjj. The occurrence of bicontinuous phases will be discussed in a later paper. 

One sees that the phase behavior is reasonable and in accord with packing considerations; as the headgroup 
increases in volume, the system passes through a series of phases from the inverted ones with most curvature, through 
the lamellar phase, to the normal ones of most curvature. We also note that as the temperature increases, the lamellar 
phase becomes unstable to one or the other of the hexagonal phases. This is very clear for the Hj phase, less so 
for the Hji phase, but is true for the latter over much of the temperature range as the slope of the La/ Hn phase 
boundary is slightly positive below T* of about 0.045. This transition with temperature is well known in the analogous 
case of diblock copolymers (Leibler, 1980). In order to model a lipid which, like phosphatidylserine, adopts the Hu 
configuration when essentially neutral (CuUis et al. 1985, Bezrukov et al., 1999), we have chosen 1/(1 -I- 74) — 0.24 in 
our subsequent studies. For comparison, the value appropriate for DOPE, calculated from the molecular volumes in 
the literature (Rand and Fuller, 1994) is 0.254. For the convenience of the reader interested in carrying out similar 
calculations, we present in Table I the values of the first three non-trivial Fourier components of the headgroup density 
0h(r), the lattice parameter D/Rg, and the free energy VtmfVh/kTV for the La, Hu and hccjj phases at T* = 0.04. 
We note, in passing, that the relative intensities of X-ray Bragg peaks can be determined directly from the Fourier 
components of the various densities with which they are associated. 

The effect on this neutral lipid of the addition of a solvent of small volume, characterized by 7^ = Vg/vh = 0.1, close 
to the value of 0.096 (Rand and Fuller, 1994; Kozlov et al., 1994) appropriate to water and a phosphatidylethanolamine 
headgroup, is shown in Fig. 2. There is a lamellar phase at small solvent volume fractions and low temperatures. 
This phase becomes unstable with respect to the Hu phase which envelops it at higher temperatures. There is a large 
region of two-phase coexistence between the ordered lipid-rich phases and an almost pure solvent phase. These features 
are reasonable, and are observed in the systems of aqueous dialkyl didodccylphosphatidylcthanolamine and of diacyl 
diarachinoylphosphatidylethanolamine (Seddon et al, 1984). Of particular interest is that we find a small temperature 
region of re-entrant hexagonal-lamellar-hexagonal transitions, an unusual feature which has been observed in DOPE 
(Gawrisch et al. 1992, Kozlov et al. 1994). As a consequence, there is an azeotrope at which the transition between 
lamellar and hexagonal phases occurs without a change in the concentration of water. We have used the coordinates 
of this point, T* = 0.06 and (f>s = 1.42, denoted Tq and {(j)s)o, to normalize the temperature and solvent-density axes. 
There is a small region of bccji in our phase diagram. Again, the possible occurrence of bicontinuous phases will be 
examined in a later publication. The uncertainty in the temperature of the phase boundaries, ST/T introduced by 
the truncation of the number of basis functions, is approximately 2 x 10^'^. 

As the volume fraction of water is increased, we find that the period of all structures increases, as is expected. In 
Fig. 3, we compare experimental results on DOPE taken in the inverted hexagonal phase at a temperature T — 22°C, 
just above that of the azeotrope (Rand and Fuller, 1994), to our values calculated just above the azeotrope. Knowing 
the molecular weight of DOPE, we convert the volume fraction of solvent, 7s<^s, which occurs in the calculation, to 
the experimental variable of weight fraction of water. The lattice parameter of the hexagonal phase in the calculation, 
however, is measured in units of the radius of gyration of either lipid tail. What value should be taken to model DOPE 
is unknown. Hence we have used in the comparison the lattice parameter D, in units of Dq, the lattice parameter at 
the azeotrope. The agreement is excellent. 

An effective value of the radius of gyration can be defined as that value which brings agreement between the 
calculated and measured lattice parameters. As the former, at the azeotrope, is D{To) = Dq — 4.79i?g, and the 
latter is 58. 9A (Rand and Fuller, 1994), the equivalent radius of gyration at the temperature of the azeotrope, Rg, is 
12. 3A for a single tail, not unreasonable when compared to the extended length of a single chain of DOPE which is 



approximately 26A. 

As the temperature of the system is lowered, the period of all structures increases, which is due to the lengthening 
of the tails as their entropy decreases. We would again like to compare our results with those of the DOPE system. 
In order to do so, we must address the temperature dependence of the radius of gyration, Rg, which appears in the 
theory, and which is not given a priori. Because the model chain is flexible, the radius of gyration is related to the 
mean square end-to-end distance, R, by a numerical constant, Rg ~ R/^/6, so their dependence on temperature is the 
same. To compare our results to DOPE, we shall assume that the temperature dependence of the radius of gyration 
which appears in the calculation is the same as that given for lipid chains by the Rotational Isomeric States Model, a 
model which describes the properties of such chains very well (Flory, 1969; Mattice and Suter, 1994). Thus wc assume 

1- <COS(h>^ ^/^ 



Rg{T)^c[—-^ T^ ' (70) 

'^ \1+ < cos (/) > J 

where the angle (j) takes the values 180° and ±70° corresponding to trans, gauche^ and gauche~ configurations, and 
c is a constant. The statistical average of cos(j> is 

cos(180°) + CT cos(70°) + cr cos(-70°) 

< cos cj) >= ^ '- -^—^ ^ '- (71) 
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with a = cxTp{—Trism /T) and Trism ~ 280.25°K. From the behavior with temperature oi Rg{T), the lattice parameter, 
D{T), at any temperature can be obtained from 

D{T) ^ [D{T)/Rg(T)] Rg{T) 

D{n) [DiTo)/Rg{To)]Rg{To)- ^' ^ 

Again, it is the factor D{T) / Rg{T) which occurs naturally in the calculation. 

A comparison of the experimentally measured (Tate and Gruner, 1989) and theoretically calculated lattice param- 
eters vs. temperature is shown in Fig. 4. In part (a), the variation of the parameter of the Hu is shown at two 
different lipid weight fractions. The agreement is very good. In part (b), the comparison is made of the Hu and La 
parameters along the coexistence with excess water. Note that this comparison is a much more stringent test, for 
it requires not only that the dependence of the lattice parameters on solvent concentration and on temperature be 
reproduced well by the calculation, but also that the phase boundaries be given well. Considering these requirements, 
the agreement is rather good. It should be noted that the agreement in 4(b)does not depend on the exact temperature 
of the triple point, which is difficult to locate precisely, but only on the existence of stable Hu and La phases which 
coexist with excess water. 

As seen in Fig. 4(b), the lattice parameter of the Hu phase is much larger than that of the La at the triple point. 
This is due to the coexistence with excess solvent, which swells the hexagonal cores, but which is only weakly present 
between the lamellae. In contrast, when the lamellar and hexagonal phases are only in two-phase coexistence with 
one another and there is no excess solvent, the hexagonal phase in general has a smaller lattice parameter than the 
coexisting lamellar phase, as shown in Fig. 5. This is due to the fact that over almost all of their coexistence region, 
the Hji phase has a smaller volume fraction of solvent than does the La phase, as can be seen from the phase diagram 
of Fig. 2. Only over the re-entrant region, which occurs as the triple point is approached, does this balance shift. 
This shift in the relative size of the latter parameters through the re-entrant region is in agreement with experiment 
(Kozlovet al., 1994). 

It is of interest to determine if any one effect can be said to drive the transition from the Hu to La phase in 
the neutral system. To this end we examine the individual terms in the thermodynamic potentials per unit volume 
flmfVfi/kTV = Evh/kTV — SiVh/kV — {SsV^/kV -\- ipg In Zs) of the La and Hji phases as the transition is crossed by 
increasing the solvent fugacity, Zg, at constant temperature T /Tq = 0.67. In Table II we show the contributions to the 
free energy per unit volume of the La phase and that of the Hu phase coming from the interaction energy, Evh/kTV, 
the lipid tails, —SiVh/kV , and the solvent —SsVh/kV — 4>s In Zg. All contributions are evaluated at the transition itself, 
which occurs at Zg ~ 3.14, and are measured with respect to the free energy per unit volume of the disordered phase. 
We also show the difference in the contribution of each term to the free energies of each phase, and the derivative 
of this difference with respect to the solvent fugacity. The difference in the contribution of the entropy of the lipid 
tails is positive because, with the lipid architecture we have chosen, the large tail volume relative to that of the head 
favors the hexagonal phase. The interaction energy favors the lamellar phase, as does the solvent, presumably because 
the interstices of that phase are two-dimensional, while those of the inverted hexagonal phase are one-dimensional. 
The difference between the lipid entropy contributions decreases with increasing solvent concentration because the 
packing constraints in the Hu phase become more severe as the size of the cores increases (Gruner, 1989). However 
it is apparent that neither this term, nor either of the others, changes so rapidly with solvent concentration compared 
to the others that any particular effect can be said to "drive" this transition. 
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VI. RESULTS: THE CHARGED LIPID 
A. Coulomb Interactions Only 

We now turn on the negative charge of the headgroups by varying the chemical potential of the free counter ions 
while enforcing charge neutrality. Increasing the density of free, positive, counter ions in our closed system is equivalent 
to increasing the magnitude of the negative charge density on the headgroups. It therefore corresponds to an increase 
in the pH of an experimental system. The charge on the headgroups is annealed, meaning that it is determined by the 
local value of the electrostatic potential, and therefore the headgroup charge varies with the location of that group. 
The parameter /?* , defined in Eq. 0, measures the strength of the Coulomb interaction. It can be written as the ratio 
of two lengths, /?* = £,/Li, where ^ = e^ jtkT is the Bjerrum length, and L; = vh/^t^I^ is a length characterizing 
the architecture of the lipid. It is reasonable that (3* be larger than or of order unity, and we have arbitrarily taken 
(3* = 1. 

We find for this system that the effect of increasing headgroup charge is to decrease the temperature of all transitions. 
As an example, we show in Fig. 6 the way the temperature, T* , of the transition between La and Hn phases varies with 
the magnitude of the average charge density of the headgroups ph = Ph;i- The range of charge density corresponds to 
the headgroups varying from being neutral to fully charged. The region beyond the almost vertical line at \ph\ ~ 0.24 
would correspond to the headgroups being charged with a probability greater than unity, and is therefore unphysical. 
The maximum value of \ph\ is slightly different in the two phases. Over the physical range, one sees a 10 % decrease 
in the transition temperature, and therefore, a stabilization of the Hjj phase with respect to the Lq, phase. Although 
at first surprising, this is a reasonable result, and can be understood as follows. In the neutral system, the single 
interaction parameter in the system, Xj is a measure of the tendency for the system to order, and is proportional 
to the inverse temperature. When the charges are turned on, the Coulomb repulsion between headgroups opposes 
ordering, and therefore has an effect similar to a decrease in x or an increase in temperature. As we have noted earlier, 
increasing the temperature causes the lamellar phase to become unstable to the hexagonal phase. A complimentary 
view is to note that increasing the Coulomb repulsion between headgroups increases the area per headgroup, and 
thus the area between headgroups, an effect which, again, is equivalent to increasing the temperature (Seddon and 
Templer. 1995). The increase in area permits the tails to splay further out, favoring the formation of the inverted 
hexagonal phase. The very same result of increasing the density of charges in a system is seen in the lyotropic phases 
of polyelectrolytes. There one finds that the phase diagram depends only on the ratio of x to the charge on the 
polymer, at least in the absence of fluctuations (Nyrkova et al., 1994). Hence an increase in the charge is equivalent 
to a decrease in x or an increase in temperature, leading to an instability of the lamellar phase. 

We are confident, therefore, that the transition from L^ to Hu phases with increasing charge on the headgroups is 
the correct result of the model. However it is completely opposite to the experimental results on phosphatidylserine 
in water (Hope and CuUis, 1980, Bezrukov et al., 1999). Hence there must be a crucial physical mechanism which is 
not included in the model which includes Coulomb interactions only between objects with a net charge. One such 
mechanism is the short-ranged attraction between the charges of the headgroup and the dipoles of the neutral water 
molecules which will cause the latter to aggregate around the former. Furthermore the counterions attracted to the 
head group will now also be associated with waters of hydration endowing them with a non-negligible volume. The 
additional volume of all these waters will increase the effective volume of the headgroup, and therefore tend to stabilize 
the lamellar phase. This additional volume will increase with the charge of the headgroup, and could be sufficient to 
counteract the Coulomb repulsion tending to stabilize the hexagonal phase. This competition is investigated next. 

B. Coulonib and Short- Ranged Solvent Interactions 

We now include in the system an attractive contact interaction between charges and the neutral solvent. This 
models the short-ranged, thermally averaged, interaction between charges and the water dipole, one which varies 
with separation r as (fcr/6)(u/e^)^(^/r)^ = uj{r), with u the dipole moment of water and ^ the Bjerrum length. 
The above expression is valid for distances such that r/^ > (u/e^)^/^. For water, ^ « 7 A, and uj{r)/kT « 4.6 x 
10^** (^/r)^. To approximate this short-ranged interaction by a contact interaction of dimensionless strength A is 
equivalent to employing Lo/kT evaluated at some fixed distance. Any reasonable choice shows that A is small. We 
have arbitrarily chosen A = 0.1, which corresponds to Lo{r)/kT evaluated at 1.9A, a distance within the regime in 
which the approximate expression for the charge-dipole interaction is valid. 

The results for the system with the Coulomb interaction of strength /3* = 1 and which includes the interaction 
between charges and solvent of strength A = 0.1 are qualitatively different from those of the charged lipids without 
this interaction. The temperatures of all transitions between ordered phases now increase with increasing headgroup 
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charge, whereas previously they decreased. Shown in Fig. 7 is the temperature of the transition between Hn and La 
phases as a function of the magnitude of the average charge density on the headgroups. One sees that with increasing 
headgroup charge, the inverted hexagonal phase becomes unstable with respect to the lamellar phase, just as in the 
experiments on phosphatidylserine in water (Hope and CuUis, 1980, Bezrukov et al., 1999) 

In Fig. 8a we show the volume fraction profiles in the La phase at a value of Zc = 0.1228, Zg = 3, and temperature 
T* = 0.04 at which the Hn — La transition occurs. The position through the system is divided by the lattice 
parameter of the phase which, in units of the radius of gyration of the lipid, is Dl/ Rg ~ 4.14. All looks reasonable. 
In particular, the volume fraction of solvent within the bilayer is negligible. This is true also when the lipid is neutral. 
In Fig. 8b, the charge density profile of the same structure is shown. One sees that the charge on the headgroup 
mimics, but does not reproduce, the headgroup volume fraction. This is because the charge on the headgroup is not 
fixed, but varies with the local electrostatic potential. The counter ion density is fairly uniform because we have 
employed a single dielectric constant, that of water, throughout the system. Were we to require that the dielectric 
constant be position-dependent, and to take the much lower value in the tail region appropriate to it, the counter ion 
density there would be much reduced. In Figs. 9a and 9b we show the same quantities for the Hu phase at the same 
value of Zc- The cut through the system is taken along the nearest-neighbor direction, and the distance is normalized 
to its lattice parameter Dh ~ 3.83i?g. Again the wavelength of the Hu phase is smaller than that of the La phase 
because, in two-phase coexistence, i.e. in the absence of a reservoir of excess water, the cores of the cylinders are 
not swollen with water and the hexagonal phase contains a smaller volume fraction of water than does the lamellar 
phase. One can infer from Fig. 9a that in the nearest neighbor direction, there is far more interdigitation of lipid 
tails than in the lamellar phase. This makes sense, as the tails must certainly stretch to fill the space between cores 
in the second-neighbor direction, so that interdigitation is expected in the nearest-neighbor direction. 

To investigate this transition further, we show in Table III the contributions in the La phase and in the Hu phase of 
the various terms in the thermodynamic potential per unit volume flmfVh/kTV = {Ei-\- E2 + E^)vh / kTV — SiVh / kV — 
{SsVh/kV + 4>s InZs) — {ScVh/kV + pc^nzc). Here Ei is the hydrophilic-hydrophobic interaction proportional to x^i 
E2 is the electrostatic interaction proportional to f3* , and E3 is the charge-solvent interaction proportional to A. 
These contributions are evaluated at the transition itself, and are measured from the free energy per unit volume of 
the disordered phase. We also show the difference between these contributions to each phase, and the derivatives of 
each of these differences with respect to the counter ion fugacity, Zc. There are several interesting things to note. 
The electrostatic energy is a relatively small contribution to the free energy of each phase, and hardly differs between 
them. Therefore it does not have a large effect in bringing about the transition. The contribution of the counter ions 
to the free energy of each phase is of the same order of magnitude as the electrostatic interaction and, like it, does not 
change rapidly with the counter ion fugacity. The contribution of the short-range charge- solvent interaction is small, 
but it changes most rapidly with the counter ion fugacity, and therefore appears to be most important in actually 
bringing about the transition itself. 

The physical mechanism in the experimental systems appears now to be clear. The lipid with an almost neutral 
headgroup forms the Hjj phase because the volume of the headgroup is relatively small compared to that of the entire 
lipid. As the charge on the headgroup is turned on, it attracts an increasing volume of waters of hydration via the 
attractive interaction between the charge and the dipoles of water. In addition, more counter ions, enlarged by their 
own waters of hydration, are attracted to the headgroup. Thus the headgroup becomes effectively larger, and drives 
the transition to the La phase. In doing so, it overcomes the effect of the Coulomb repulsion between headgroups 
which, in fact, opposes the transition to the La phase. 

We have argued earlier that an increase in Coulomb repulsion was similar to an increase in temperature, and 
therefore its effect in favoring the Hjj phase could be understood from the phase diagram of Fig. 1 . In a similar way, 
an increase in the strength of the charge- solvent interaction increases the volume fraction of the headgroup, so that 
one moves to the right in Fig. 1 and stabilizes the lamellar phase. The combined effect of increasing the headgroup 
charge is to increase the effective volume fraction of the headgroup and effective temperature. The phase which is 
ultimately stabilized results from the competition of these two effects. By changing the strength of the Coulomb 
interaction, we have verified that one can alter the effect of this competition, and not only stabilize the inverted 
hexagonal phase, but also bring about the inverted b.c.c. phase or even the disordered phase, scenarios which are 
again understandable from our discussion and the phase diagram of Fig. 1 . But that the result of this competition 
can be the stabilization of the lamellar phase at the expense of the inverted hexagonal phase has been demonstrated 
by our model calculation, as well as by experiment. 
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Table I Anhydrous, neutral lipid: the lattice parameter, the free energy, and the first 3 non-trivial Fourier components 
of (/)/i(r) for La, Hu and bccu phases at T* = 0.04 and 1/(1 + 74) = 0.24. Recall that 0/i,i = (j)t,i = 0.24 for all 
phases. 





D/R, 


Vhn,nf/kTV Cf,h.2 


Cf>h,3 X 102 


^hA X 10^ 


ia 


2.921 


0.7969 0.2272 


6.781 


-1.171 


Hii 


3.167 


0.7936 0.2368 


0.476 


-2.372 


bccii 


3.400 


0.8027 0.2204 


1.094 


-4.038 



Table II Neutral lipid: contributions to the free energy per unit volume and temperature in the La phase and in 
the Hji phase, the difference in these contributions, and the derivative of this difference with respect to the 
solvent fugacity. All contributions are evaluated at the La, Hu transition occurring on the water-poor side of 
the azeotrope. The solvent chemical potential at the transition is Zs ~ 3.15. All contributions are measured 
with respect to the free energy of the disordered phase, (D). The temperature, T/Tq = 0.67. £ = VhE/VkT, 
Si = VhSi/Vk, and Ss = VhSs/Vk + (ps In Zg. 





La-D 


Hii -D La- Hu d{La - 


- Hii)/dzs 
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-0.6923 


-0.6086 -0.0837 


0.0244 


-Si 


1.1367 


0.9123 0.2250 


-0.0413 


-Ss 


-0.6666 


-0.5253 -0.1413 


-0.0394 



Table III Charged lipid with both Coulomb and charge-solvent interactions: contributions to the free energy per 
unit volume and temperature in the La phase and in the Hu phase, the difference in these contributions, 
and the derivative of this difference with respect to the counter ion fugacity. All contributions are evaluated 
at the La, Hu transition, z^ = 0.1228, and are measured with respect to the free energy of the disordered 
phase, (£>). The temperature, T* ~ 0.04, Zs = 3, £i = VhEi/VkT, with i = 1,2,3 being the hydrophilic- 
hydrophobic. Coulomb, and charge-neutral interactions respectively. Si = VhSi/Vk, Sg = VhSs/Vk + (ps In Zs, 
and Sc = VhSc/Vk + pc In Zc, 

La — D Hu — D La — Hu d(La — Hii)/dzc 



£1 


-0.6904 


-0.6059 


-0.0845 


0.0152 


£2 


0.0037 


0.0024 


0.0013 


0.0196 


£z 


-0.0378 


-0.0268 


-0.0110 


-0.1082 


-Si 


1.1373 


0.9110 


0.2263 


-0.0257 


-Ss 


-0.6365 


-0.5034 


-0.1331 


0.0215 


-Sc 


0.0030 


0.0020 


0.0010 


0.0153 
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Figure 1 Phase diagram of the neutral hpid as a function of dimensionless temperature, T* = l/2xN, and relative 
headgroup volume, 1/(1 + 74). In addition to the disordered phase, D, there are normal and inverted body- 
centered cubic phases, bccj and been, normal and inverted hexagonal phases. Hi and Hu, and the lamellar 
phase La- 
Figure 2 Phase diagram of a neutral lipid with 1/(1 + 74) ~ 0.24 in a solvent with 7s = 0.1 as a function of 
temperature, T/Tq, and fraction of solvent, 4>sl {4's)o, with Tq and ^sift's)!:) the temperature and volume fraction 
of solvent at the azeotrope. 

Figure 3 Comparison of theoretically calculated and experimentally measured values of the lattice parameter D/Dq 
of the Hi I phase at a temperature just above the azeotrope vs. weight fraction of water, (f>^. The lattice 
parameter at the azeotrope is denoted Dq- 

Figure 4 Comparison of theoretically calculated and experimentally measured values of the lattice parameter, 
D{T)/Da, vs. absolute temperature T (a) for two different weight fractions of lipid, 4>T ' ^.nd (b) along co- 
existence with excess water. The absolute temperature of the azeotrope is denoted Tq. 

Figure 5 Lattice parameter of the Hn and La phases along their mutual coexistence as a function of absolute 
temperature. The absolute temperature of the azeotrope is denoted Tq. 

Figure 6 Transition temperature between La and Hn phases for the same lipid in Figure 2, but now with a charged 
headgroup. The dimensionless strength of the Coulomb interaction is /?* = l.Thc transition temperature is 
plotted as a function of the magnitude of the average charge density on the hcadgroups. The solvent fugacity is 
fixed at Zs = 3. 

Figure 7 Transition temperature between La and Hn phases for the same lipid as in Figure 6, but now including 
the interaction between charges and neutral solvent of strength A = 0.1. The solvent fugacity is again fixed at 
Zs = 3. 

Figure 8 a) Volume fraction distribution in the La phase of the solvent, headgroups, and tails, in the above system 
at a counter ion chemical potential of Zc — 0.1228 corresponding to the La — Hn transition. The temperature 
is T* = 0.04, and the lattice parameter of the lamellar phase is Di/Rg — 4.14; b) charge densities arising from 
the headgroups, the counter ions, and the total charge density in the La phase under the same conditions as in 
a). 

Figure 9 a) Volume fraction distribution in the Hn phase of the solvent, headgroups, and tails, in the above system 
at a counter ion chemical potential of Zc = 0.1228 corresponding to the La — Hn transition. The temperature 
is T* — 0.04, and the lattice parameter of the inverted hexagonal phase is Dn/Rg = 3.83. b) charge densities 
arising from the headgroups, the counter ions, and the total charge density in the Hn phase under the same 
conditions as in a). 
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